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Abstract 

In dense quark matter, the response of the color superconducting gaps to 
a small variation, 5/i, in the chemical potential of the strange quark was 
studied. The approximation of three massless flavors of quarks and a gen¬ 
eral ansatz for the color flavor structure of the gap matrix was used. The 
general pole structure of the quasi-particle propagator in this ansatz is pre¬ 
sented. The gap equation was solved using both an NJL interaction model 
and perturbative single gluon exchange at moderate densities and results are 
presented for varying values of 5/1. Quantitative and qualitative differences 
in the dependence of the gaps on were found. 


1 Introduction 


The physics of strongly interacting matter at high densities and low temper¬ 
atures has been the subject of much research in recent years. It has long been 
known im] that at sufficiently high densities a system of quarks should form a 
condensate of Cooper pairs which breaks the SUc{^) symmetry and becomes 
a color superconductor. The formation of the condensates leads to gaps in 
the quasi-particle mass spectrum. The authors of |I| estimated that the gaps 
were of the order of A ~ where fi is the quark chemical potential. 

Recently it was shown at realistic values of /i in an instanton induced NJL 
model that gaps of the order of /i could be obtained|^. This stimulated a 
great deal of research^ in the ensuing years and has resulted in a proliferation 
of predicted superconducting ground states. These states may be realized in 
the cores of neutron stars and lead to observable effects^ 

It is widely accepted 10 that the color superconducting ground state 
at asymptotic densities is the Color Flavor Locked (CFL) state[^]: 

- 4 K) + y). (1) 

where the Greek indices are color indices, the Latin indices are flavor indices 
and the 3 and 6 subscripts refer to anti-triplet or sextet conhgurations in color 
and flavor spaces respectively. At lower densities (/r ~ m^/4A), it is likely 
that the ground state is a superconducting state involving the condensation 
of Cooper pairs in the u — d sector only (2SC). Finally at still lower densities 
the favored ground state will be ordinary hadronic matter. Two new phases 
have recently been predicted: Crystalline Color Superconductivity 1 and 
CFL with meson condensation[0. These predictions, while not necessarily 
at odds with one another, indicate that the transition region between the 
CFL state and hadronic matter is not completely understood. 

Including the strange quark mass in the gap equation with no approxi¬ 
mations introduces two sets of complications: 1) massive quarks means that 
there are 4 new types of Dirac structures allowed for the condensates, and 
coupling between the corresponding gaps; 2) the fact that the strange quark 
is different from the other quarks means that gaps involving the strange 
quark should be different from those with zero strangeness. 

In order to understand the implications of these two complications it is 
useful to separate them and understand them individually before tackling 

^For extensive references see the review articles Q. 
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the full problem. Therefore, in this paper, all three quarks are treated as 
massless and the strange quark is distinguished by giving it a different chem¬ 
ical potential (ps = p — h/i) than the u and d quarks. This is motivated by 
observing|T^ that for hxed Fermi energy, the Fermi momentum of a massive 
particle is less than the Fermi momentum of a massless particle: 


Ep — Pfo — po 


for massless particles, (2) 


! - 1 TTl^ 

Ep = \Jpp^ + m? K. pp^ -|- --= /io for massive particles. (3) 


2pF, 


This approximately corresponds to a shift in the Fermi momentum or the 
chemical potential: 


1 m? 

f^m PFm ~ 7) ~ /^O Ofl. 

^PFm 


( 4 ) 


This approach is similar to that taken in [|^] where this approach was studied 
using a four fermion interaction (NJL) model. Qualitatively similar results 
to those of were obtained in using an NJL interaction and a different 
prescription for inclusion of strange quark mass effects. In this paper an NJL 
model is considered in Section 3 for comparison to previous results and to 
illustrate the contrast with the results of the subsequent section. In section 
4, perturbative single gluon exchange is considered for the interaction of 
the quarks. This approach is valid in the high density regime and gives an 
approximate model in the moderate density regime which is distinct from 
and can be compared to the NJL interaction model. 

Understanding the effects of a shift in the chemical potential of one of 
three quark flavors is another motivation for this research. The papers p, 
|T^ |T^ |T^ also deal with shifts in the chemical potential in an NJL model, 
mostly in the two flavor case and find some interesting effects. The three 
flavor case was discussed in P and but was not analyzed in detail. In 
it was predicted that the number densities of the 3 quark flavors are 


13 


constrained to be equal up to J/r = Agg/v^, at which point there is a hrst 
order transition to some less symmetric phase of quark matter. It would 
be interesting to test this prediction using perturbation theory. Finally in 


151 hrst order perturbation theory is used to study the non-degenerate two 


havor case. In this work, hrst order perturbation theory is applied in the 
non-degenerate three havor case. 
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Perturbation theory applied to Color Superconductivity has been the sub¬ 
ject of much important research^ It has not previously been applied to the 
case of color superconductivity with three flavors of non-degenerate quarks. 

Recently numerical results were obtained |T^ in the 2 flavor case in per¬ 
turbation theory with single gluon exchange and massless and degenerate 
quarks. The Sfi = 0 results presented in this paper for perturbative single 
gluon exchange with three flavors are qualitatively similar and can be con¬ 
sidered as independent confirmation of as well as an extension of them 
to three flavors. 

The fact that the three flavor case is different from the two flavor case 
for non-degenerate flavors should be emphasized at this point. This is clear 
from the NJL analyses of the two different cases. In Q the gap is shown 
to be independent of 6fi up to some critical value. In [|^ |^, h/i has a 
signihcant effect on the pattern of condensation even at small values of 6fi. 
This difference is simply a result of the fact that the algebra of the generators 
of SU{2) is much simpler than the algebra of the generators of SU{3). The 
important point that this illustrates is that results from the two flavor case 
may not tell the full story in the three flavor case. 

The goals of this work are threefold. The hrst goal is to understand the 
effects of the non-degeneracy of the strange quark in a simplihed model of 
QCD as a hrst step toward understanding the effect of the strange quark 
mass in a more complete way. The second goal of this work is to understand 
the effects of a small shift in the chemical potential on solutions of the gap 
equation for comparison with results in other approaches to the problem. 
The third goal is to determine how hrst order perturbative results compare 
to results obtained using the NJL interaction model. 

Results presented in this paper are the poles of the quasiparticle propa¬ 
gator in this approach in a more general ansatz than in and in a diherent 
model than . Numerical solutions for the gaps at moderate densities using 
both an NJL interaction model and perturbative single gluon exchange are 
presented to illustrate the ehect of J/i and the qualitative and quantitative 
diherences in the results. 


^See for example ^ and references in |^. 
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2 Gap Equation 

The gaps resulting from the formation of a color superconducting ground 
state can be determined by solving the mean held gap equation[^]: 

A,B 2 2 

which is specialized to an interaction with the color structure of single gluon 
exchange. — q) is the gluon propagator, are the Gell-Mann matri¬ 

ces, G^i^q) is the bare antiquark propagator, G+(g) is the full quasiparticle 
propagator and A(fc) is the gap matrix. 

The gap matrix is a matrix in color, havor and spinor space that contains 
all the specihc gaps. Restricting consideration to spin zero gaps in the mass¬ 
less case the Dirac structure of the gap matrix can be written in a basis of 
four projectors!^: 

Ptik) = , e. A = ±1 (6) 

which are products of projectors onto positive and negative energy states and 
positive and negative helicity states. 

The color havor structure of all objects will be given as 9 x 9 matrices 
consisting of a 3 x 3 matrix in havor space whose components are 3x3 
matrices in color space: 
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The most general gap matrix can be written as: 

A(t)= E (Ay)'p;(fc). (8) 

e,/l=dil 
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The gluon propagator, -D^^, is the object that differs between the NJL 
model and perturbation theory. In the NJL model the interaction is modeled 
by a 4 fermion interaction and the gluon propagator is replaced by an identity 
matrix. This constitutes a low energy effective model of the interaction. 
In perturbative single gluon exchange is the bare gluon propagator. 

Perturbation theory is definitely valid at high energies and densities. At 
moderate densities perturbation theory can be considered a model for the 
quark interactions. 

The bare propagator for the a quark of flavor / is given by: 

(<^o(^))^ = (7''^..±/^/7o)“^ = ((fco±/^/)7o-7-^) ^ (9) 

If all fif are equal the particles are degenerate and the bare propagator is 
the identity on the color-flavor space. In the case considered here the bare 
antiquark propagator can be written in the form: 

G'o = ((^0 -- 1^1 7 • fc) (10) 

_ ^ ^ (p<"> + pf) p,;W 

^ {ko-n) + e\k\ {ko-n-6n5i3) + e\k\ 

where the PO) qj-q projectors onto a specific flavor sub-space, ie: 



/ 1 0 0 \ 
0 0 0 
V 0 0 0 / 


( 11 ) 


where bold font indicates 3x3 matrices in color space, and similarly for 


id) 


and P 


(d 


/ 


The quasiparticle (quasi-antiparticle) propagator is determined from the 
bare particle and antiparticle propagators and the gap matrix by the relation [PO 


G^ = \\G 


■y±l- 




( 12 ) 


where A"*" = A and A = 7 oA"^ 7 o. 

The ansatz for the color-flavor structure of the gap matrix used in this 
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research is: 



/ a+/i+e 
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0 
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where: 










a + h — e 
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^66 

{ud) + A 33 (Md), 




e 

= 

^66 

{ud) — 

^33 ('WC?)) 




c 

= 

^66 

{us) + A 33 (ms), 




/ 

= 

Aee 

{us) - 

^33{us), 




a — h = A66(ss) 


( 13 ) 


(14) 

(15) 

(16) 

( 17 ) 

(18) 


The form of this ansatz is essentially the same as in with minor differences 
to simplify the poles of the quasi-particle propagator. The correspondence 
between the different conventions is: 


b = 


a + h — e 


d = a — h 


(19) 


and e,c,/ unchanged. 

In the case of degenerate color-flavor locking: 


h 


c = 


e = / 


a — h = 


a + h + e 


( 20 ) 


In order to simplify the calculation, following the example of |^, the 
ansatz ([I3|) is rotated to a different basis. The columns (rows) are referred to 
by the combination of color and flavor indices, The {u,r) and {d,g) 

basis vectors are rotated to new basis vectors l/\/2 {{u,r) =f {d, g)) and all 

^Therefore for instance the 2nd column is the (z, a) = {u, g) column. 
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the other basis vectors are left unchanged^. Note that this change of basis 
leaves cind operators unchanged which means that the bare 

propagator has the same form in this new basis. In this basis the gap matrix 
now has the form: 




/ e 0 0 0 0 

0 0 0 e 0 

0 0 0 0 0 

0 e 0 0 0 

0 0 0 0 a + h 

0 0 0 0 0 

0 0/0 0 
0 0 0 0 0 

y 0 0 0 0 y2c 


0 0 0 0 \ 

0 0 0 0 

0/0 0 
0 0 0 0 

0 0 0 y2c 

0 0 / 0 
0 0 0 0 

/ 0 0 0 

0 0 0 a — h j 


( 21 ) 


Of course all other objects with color or flavor structure must be transformed 
to the same basis. These details are only included in order that the quasi¬ 
particle propagator can be presented below in a fairly concise way. 

Up until this point no assumptions about the Dirac structure of the gap 
matrix has been made. For ease of calculation in the perturbative analysis 
the gap matrix is assumed to have the usual form: 


A(fc) = {i(-4)C759f(4)> = Ay ehPi(k) (22) 

e,h=^l 


Substituting (^2]) and (|^) into (|T^, the quasiparticle propagator is: 


G+ = 


/ P 0 0 0 0 0 

0 P 0 0 0 0 

0 0 Pi 0 0 0 

0 0 0 P 0 0 


0 0 0 


0 0 0 0 0 0 

0 0 0 0 0 0 

V 0 0 0 0 P 0 


0 

0 

0 

0 

0 


0 

0 


0 

0 

0 

0 


0 A I 

0 0 0 0 0 Pi 0 

P2 


0 \ 
0 
0 
0 


0 c 
0 0 


0 

P2 

0 


0 
0 
B j 


(23) 


^This means for example that the new first column is 1 /v^ ((u, r) — (d, g)) and similarly 
for the fifth column. 
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where: 


E = 

Fi = 

F 2 = 

A =- 


Qo + I 


{Qo - 0 ('?o + 0 - 

_ qo + l -6fi _ 

(?o - 0 (<?o + ^ - - P 

Qo + ^ 


A+ + [ / ^ -/' ] A“ 

A+ + [ / ^ -/' ] A- 
A+ + [ / ^ -/' ] A- 


(go - ^ + Sjj.) (go + 0 - P. 

(a — h)‘^{l + gO) + {2p + (/ + go)(^ — Qo ~ F Qo ~ <^h))) 


(24) 

(25) 

(26) 
A+ 


(g§ - e+(a, h, c, I, (gg - e_(a, h, c, /, 5/i)2) 

- [ / ^ -r ] A- (27) 

(a + h)^(/ + gO — 5/i) + {2p + (/ — go)(/ + go ~ ^a))(^ + 5'o))) 


B = - 


C = 


{ql - e+(a, h, c, /, /r, (5/i)2) (gg - e_(a, h, c, /, /i, 5/x)2) 
[ / ^ -/' ] A- 

■\/2c (a(2/ + 2go — 5p) — h5p) 


{ql - e+(a, h, c, I, /i, 5/i)2) (gg - e_(a, h, c, /, /r, 5/i)2) 

+ [/^-riA- 


A" 


and: 


e±(a, h, c, /, 6fi) = l{l — 5p) + 5p + a + 2p + h 


A+ 

(28) 

(29) 

(30) 


± + h^)+2ah(2l - + (21 - l>,f) , 


using the dehnitions: 


I = (k|-h), (31) 

= (kl+h)) (32) 

which are the only ways that the functions depends on |g| and /i. 

The form (p3|) of the propagator is consistent with expectations. The 
strange quark has been distinguished from the up and down quarks. The 
ansatz still has flavor locked to color and therefore the blue quark has been 
distinguished from the red and green quarks. The (u, g) and (d, r) quarks 
propagators are unchanged. The (n, h) and (d, h) propagators are affected in 
exactly the same way. The (s, r) and (s, g) propagators are also affected in 
exactly the same way. If one rotated back to the original basis one would 





















find the (m, r) and (d, g) propagators are also affected in exactly the same 
way. The propagator could be completely diagonalized to separate the poles 
which would simplify the calculation in some ways but would complicate it 
in others, so this form is used in this work. 

The poles of the quasi-particle propagator are therefore: 


go = -p £2 

degeneracy = 3, 

(33) 


degeneracy = 2, 

(34) 


degeneracy = 2, 

(35) 

gO = ±e_ (a, h, c, 1, Sg) 

degeneracy = 1, 

(36) 

gO = ±e+(a, h, c, 1, 6g) 

degeneracy = 1. 

(37) 


It can be shown that in the case where the sextet gap^ are neglected, the 
poles of this propagator reduce to those given in . Further the hrst eight 
of the poles reduce to the octet poles and the last pole reduces to the singlet 
pole at d/i = 0. The poles obtained in this model are based on a different 
model than |T(] . 

Substituting the ansatz (^), the solution for (7'*~ (p3D , and the bare prop¬ 
agator into the gap equation one obtains a matrix gap equation. 


A{k) 


-ig^ j “ ?) 7 /. 7 ° 7 ^ 7°7 ■ k 

(^A^(g)M(a, h, c, I, go, d/i) — A~{q)M{a, h, c, - 

-ig^ J ■ ^7° 

{h.~{q)M{a, h, c, I, go, d/i) - A~^{q)M{a, b, au, 


(38) 

I', go, d/i)) 

(39) 

-l',qo,Sfi)) 7 ^ 


®The sextet gaps are much smaller than the anti-triplet gaps. 
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where D^^{k — q) = SAsD^'^ik — q) has been used, the index A has been 
summed over and the A matrices have been multiplied through. 


The matrix M( 

a, h. 

c, 1, qo, Sqi) has the form: 






( Mn 

0 

0 

0 

0 

0 

0 

0 

0 

\ 



0 

0 

0 

Mn 

0 

0 

0 

0 

0 




0 

0 

0 

0 

0 

0 

M 37 

0 

0 




0 

Mn 

0 

0 

0 

0 

0 

0 

0 



M = 

0 

0 

0 

0 

M 55 

0 

0 

0 

M 59 


(40) 


0 

0 

0 

0 

0 

0 

0 

M 37 

0 




0 

0 

CO 

0 

0 

0 

0 

0 

0 




0 

0 

0 

0 

0 

M 73 

0 

0 

0 




V 0 

0 

0 

0 

M 95 

0 

0 

0 

Mqq 




which should be compared to the form of the ansatz ( 0 ). In general M 73 7 ^ 
M 37 and M 95 7 ^ M 59 which can be shown to be inevitable if one examines the 
product Gg carefully. However, the terms which spoil the symmetry 

of these components can be shown to vanish under the qo integration over the 
interval (—cx),cx)). This is immediately clear in the NJL model and follows 
with a little more care in perturbation theory. Therefore the matrix gap 
equation does close under this ansatz. 

Everything discussed so far follows without specihcation of the interac¬ 
tion model. In the following section the gap equations in the NJL interaction 
model will be given for the ansatze above. Solutions for the gaps will be pre¬ 
sented in order to compare to previous results and to highlight the different 
behavior of the gaps using perturbation theory. 


3 NJL Interaction Model 

In an NJL model the fermion interaction is taken to be: 

^ Gg>^’'6AB 


(41) 


which is the same form as the interaction model of and the gap equation 
becomes: 


A{k) = -iG 


d}q 

W)‘ 


7 ° 7^7 ■ k'^^ 


(42) 


(^A {q)M (a, h, c, /, gg, J/i) - A+(g)M(a, h, c, go-<5/i)) 'Ju 
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Using the relation: 7 '^ 7 ° 7 * 7 ;y = 0, multiplying by 7 ^ and tracing over 
Dirac indices one obtains: 

A“/ = AiG j {M{a, h, c, /, go, + M{a, h, c, -I', go, 5p)) (43) 

where the right hand side of the equation does not depend on k anymore and 
the functional dependence on the left hand side has been dropped. Equating 
coefficients of matrices on each side produces a set of coupled gap equations: 


'-“Iwpf- 


+ [I ^ -/'] 


+ [/ ^ -/'] 


a = AiG 


2(a - h)c^ - (g + h) ((a - hf + (/ - - gg) l , 

4(g2-e^)(g2-e2) J 

^ A-r^ f A [c(a2 - 2 c 2 - -/2 + g02 +/5/i) 

(27r)4[ 2{gi-el){qi-el) 

/ / 1 / 

12(/2 + (/ + g0)(/ - gO - (5/i)) 12(/2 + (/ - g0)(/ + gO - Sfi))_ ^ ^ 

“ = / (liy 78(e2 + P 

2(5a + 3h)c^ + (a + h){2l — 5jj,)5iJ, — (5a — 3h){P — g^) — {a — h){a + h){5a + 3h) 


24 {ql - el) (gg - e^.) 

+ [/ - -/'] 

(47) 

771 _ 8 (e 2 +% - qi) + 

2(3a + 5h)c^ — (a + h){2l — 5iJ,)5jj, — (3a — hh){P — ql) — [a — h){a + h){3a + 5h) 


24(gg-e^) (gg - ) 

+ [/ ^ -/'] 

f d^q [c(a2 - 2 c 2 - -/2 + g02 +/5p) 

' = J (2^7 [-6(7-4)M-eO- 

/ , / 1 , r, . 


(27r)4 

/ 


4(/2 + (/ + g0)(/ - gO - <5p)) 4(/2 + (/ - g0)(/ + gO - 6^i)) 


+ [I ^ -/'] 


These gap equations include the terms I —I' which correspond to the 
contribution of the anti-particle gap and do not lead to a gap at the Fermi 


11 


surface. The anti-particle gap is not necessarily small and in fact is equal 
to the particle gap in the NJL model. These terms are suppressed, however, 
by at least 1/p^ and are usually neglected as is done here. This is crucial 


because the antiparticle gap is gauge dependent (see for example |TT[| ). 

Evaluation of the integrals on the right hand side of the gap equations 
is facilitated by the analytic continuation, go The ^4 integral is 

done by contour integration closing the contour in the upper half plane and 
picking up the poles, (^ ) -(|37| ) , which lie in the upper half plane after analytic 
continuation. The angular integrals can be done trivially reducing the right 
hand side of the gap equation to an integration over g(= |^) which can be 
done numerically. This approach has the advantage that the quasi-particles 
are automatically on the mass shell. 

The range of integration for q is not inhnite since the NJL model is a 
four-fermion interaction model and must have an UV cutoff. In this work 
the integrals are simply regularized by the factor: 


n(q) = 


h* 


(q^ + A^y- 


( 60 ) 


with A = 1000 MeV chosen so that the magnitude of the results in the NJL 
model are comparable to the perturbative results. 

The coupled gap equations were solved iteratively for p = 500 MeV using 
numerical integration and the results are shown in Figure |l|, after conversion 
to the convention of for comparison. The effect of a shift in the strange 
quark chemical potential is a disruption to the pairing of the heavy quark with 
the light quarks in the dominant channel and a slightly smaller enhancement 
of the pairing of the light quarks. The results shown in Fig. exhibit the 
same qualitative behavior as shown in Figure 6 of [|T^]. Note that the range 


of J/i shown in Fig. does not go all the way up to the phase transition 
illustrated in [|^ and the shifts seen in the results of this paper are larger 
than in 
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The results are shown in terms of gaps with specihc color-flavor symmetry 
in Fig. 1^. These results show that the Agg gaps are non-zero but small. As 
well, the effect described above is seen more clearly: significant disruption to 
the pairing of the heavy quark with the light quarks in the dominant channel 
and a slightly smaller enhancement of the pairing of the light quarks. These 


results qualitatively agree with the results of |]^ as far as the increase in 
A^^{ud) and decrease of in A 33 (ms). The quantitative agreement appears 
at hrst sight to be rather worse than above. The differences between the 
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5 |X [MeV] 


Figure 1: Numerical solutions for gaps as a function of 5/1 in the conventions 
of 0 . 


results of [O and |]OT, however, are mainly the result of different parame¬ 


ter choices pin . Therefore the results presented here for the NJL model are 
consistent with previous analyses. 

The mathematical reason for this phenomenon can be seen by examin¬ 
ing the gap equations for the individual A 33 gaps, but neglecting the Age 
gaps and the difference between the Agg gaps on the right hand side of the 
gap equation. Expanding the integrands on the right hand side of the gap 
equations one obtains: 


^33 M = 4 :G 
^ l{2f 


+ 5Agg 


dq 
;) 


n{q) 


2 A 33 


^33 


%(P + + 9(!^ + 4Ay 

H2e + i4Ay \ 


162Agg(/2 + A|)3/2 162Agg(/2 + 4A|)3/2 


' 8 /® + 28/^AA + 29/2 AF + 18A| 


+0{6^/‘ 


^33 -r ^33 -r 

486A|5(P + Ai3)V2 


33 


8/® + 88/^ AA + 314/2 AF ^ 395^6 ' 

l86AyFT4Ay575 , 


(51) 


5/i2 
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6 |X [MeV] 


Figure 2: Numerical solutions for gaps as a function of djj, with specific color 
flavor symmetries. 


A 33 (ms) =AGj^q^ lZ{q) 

( 1{P + TAgg) 


2A33 A33 

9(/2 + AyV2 +9(/2 + 4Aig)V2 

/(/2 - 2 Ay /A 33 


162A33(/2+Ay3/2 162A33(/2 + 4Ay3/2 12(/2 + Ay3/2 

2/® + /^A|g + 8/2A|g - 9A|g 2/® + 16/^A|g + 14/2A|g _ 

486AyFTAyv5 486Ay/= + 4A|,)V2 


V 486A|3(P + A|3)V2 486A|(P + 4A|3)V2 

, (2/^ — Agg)Agg^ ^2 , /o/r 3^ 

^48(/2 + Ay®/2j'^^ J ^ 

One expects these equations to be dominated by the region near the Fermi 
surface, I = 0 {q = fi). Expanding about / = 0 one obtains an approximation 
for the integrand in this region of the integration: 


A33(-Uc/) : g 2 7 ^(g) +C>(( 5 y/Agg)+C>(//A33) ( 53 ) 


648 AA 


432 AA 


^ssius) : g2 7^(g) 


83 / (5/i 


27 AA 


+ 0(5y/AA) + 0(//A33). (54) 
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The I > 0 part of this approximation is dominant because of the factor, 
= (/ + /i)^. Therefore, if this region of the integration is dominant, A^^{ud) 
should increase linearly with 5/i and a quadratic component for larger values 
of 6fi. Similarly A^^{us) should decrease linearly with 5/1 and a quadratic 
component for larger values of 6fi. The decrease in A^^{us) should be greater 
than the increase in A^^{ud). This is exactly the type of behavior exhibited 
by the solutions as shown in Fig. 

As was mentioned above, this behavior was already seen in a different 
model in |^. One reason for presenting these results is as a check of the 


method to show that the model used here gives results that are consistent 
with previous analyses. The second and more important reason is to contrast 
with the results in the next section using perturbation theory and to explain 
the differences. 


4 Perturbative Single Gluon Exchange 


In the perturbative analysis, following the analysis of [0, the gluon propa¬ 
gator: 


Df,uiq) = 


P. 


_ 


+ 


P. 


fll' 


-Ul ' - Hriq) 

is used in the weak coupling limit where: 




qi^qi^ 




A A ^^^.7 n n 
q 4 ~ QiQj' 


and using: 


with: 




■2 

rriD = 


9 = 4ir 


2 ir 2 

12ir 


(lliV,-2iVy)log 9 VA^cd 


(55) 

(56) 

(67) 

(58) 

(59) 


Evaluating the right hand side of the gap equations is more complicated 
in this case. The gluon propagator depends on the angle, 6, between the mo¬ 
menta of the scattered fermions and this integral must be done numerically. 
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There is still one trivial angular integral that can be done immediately. As 
in the previous section, the analytic continuation go Ia is done and the 
Qi integration is done by picking up the poles (|33|)-(|T7D of the quasiparticle 
propagator in the upper half plane. The right hand side of the gap equa¬ 
tions is then reduced to integration over g(= |^) and one angular integral, 
/ dcos6, which can be done numerically. This approach has the advantage 
that the quasi-particles are automatically on the mass shell. 

The gluon propagator is not the same for each term since it must be 
evaluated at the poles of the quasiparticle propagator. The poles of the 
gluon propagator are ignored in this calculation. 

As was done in the previous section and in |T^ the contribution of the 
antiparticle gaps is ignored. This is because the antiparticle gap does not 
lead to a gap at the Fermi surface and it’s contribution in the gap equation 
is suppressed by at least a factor of As well, in a recent paper 


m it 


was shown that the anti-particle gap is further suppressed by an extra power 
of the coupling in perturbation theory if the quasiparticles are on the mass 

shell. 

In order to make a numerical solution of this problem feasible we assume 
that the dominant contribution to the gap equation comes from the mean 
value of the ud and us gaps which are antisymmetric in color and flavor: 


A A 33 (ud)-h A 33 (-us) 1 1 

^33 =-^-= g(a + ^-3e) +-(c-/). (60) 


The other gaps and the splitting between A^^{ud) and A^^{us) were evaluated 
after solving to check that they were significantly smaller than A 33 and could 
safely be neglected. 

The gap equation in this case is: 

A33(p) = dq j_^dx ^f^gf F(g,A 33 (g),d/i) (61) 

3/2-1/2X l/2 + l/2a; 

p 2 _|_ q ,2 _ 2 pqx -\- [If-r — (p4 — g4)^] — 2pqx + \n.L — {Pa — ^ 4 )^] 

where g* are the poles of the quasiparticle propagator, (|3^ - (|37D and x = cos 6. 
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(62) 


Fiq,^33{q),Sn) = 271 




+ 


-^aa 


P + AF 12^{l-Sf7/2y + AF 


^aa ((2^ + Six)5ix + 7 A 33 + D(q, A^Aq)^ 

24:y/2D{q, A^Aq), Sp) ( 2 P - 2W/i + SjA + SA^g + D{q, A^Aq), Sp)) 
^33 ((2/ + 5/i)5/i + 7Ag3 — D{q, AsAq)^ ^p)') 


2A^/2D{q, AiAq), ^A) - 2W/i + 5A + 5AF - D{q, A^Aq), 


D{q, A-sAq), 5A = ^9A|3 + 2AlApiAl + 3(5/i) + 5AiAl - 6A (63) 

One can solve the gap equation (^Tj) iteratively using numerical integra¬ 
tion with some care. In particular the collinear singularity at \p\ = \q\ and 
a; = 1 in the gluon propagator must be handled. In addition there are other 
features of the integrand near x = 1 which are nonsingular but must be 
handled with a little care. 

Perturbative analysis is certainly reliable at asymptotic densities where 
the Fermi momentum is very high. If color superconducting quark matter 
is actually observed, however, it will be at moderate densities of 4-5 times 
nuclear matter density in the core of neutron stars. At /i = 500 MeV where 
the results below are obtained, the perturbative single gluon exchange in¬ 
teraction should be considered a model which is distinct from and can be 
compared to the NJL model which is a low energy model. 

Using this procedure a solution of (^) for A^Ap) obtained for differ¬ 
ent values of 6 fi. These results were then used to calculate estimates for the 
individual gaps. The results are shown for 6 fi = 0 MeV and = 45 MeV 
in Figure The error in the solution for A 33 is estimated to be less than 
0 . 01 % and the error in the values for the individual gaps are probably less 
than 0.25%. 

Comparing the = 0 results to previous analyses one sees that the 
peak value of A^Ap) is of the same order of magnitude as predicted by the 
analysis of where A 33 (po) is determined for the case when both magnetic 
and electric gluon exchanges are taken into account. The form of the Sp = 0 


results is qualitatively similar to results found in [|T^ for A^Ap) i^ fh® two 
flavor case at a similar density. 
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Figure 3: Numerical solutions for the gaps as a function of p for = 0 MeV 
and (5/i = 45 MeV. 


The agreement with |]^ includes a small cusp feature seen most clearly in 
the (5/i = 0 results in Figure ||. This feature arises because the magnetic glnon 
propagator is non-analytic at po = go dne to the form of nr(p —g) bo-goI 
in (|57D . This leads to a non-analyticity of the gap fnnction as a fnnction 
of p. This featnre arises in the results of |]^ as seen in their Fig. 6 (a). It 
is a small effect, bnt the fact that it occurs in both analyses snpports the 
conclusion that the independent analyses are consistent. 

The neglected Age gaps are less than 5% of the A 33 gaps, and A^^^ud) 
and A 33 (ms) are within 5% of their mean valne, A 33 , at the maximnm valnes 
of 6fi investigated. This indicates that the estimates obtained by solving for 
A 33 and then calculating the individnal gaps shonld be fairly reliable. 

The most important result of this analysis is that both A^^{ud) and 
A^^{us) decrease with 6p. This is in contrast to the resnlts obtained in 
the NJL model where only the A^^{us) gap decreases and the A^^{ud) gap 
increases. The peak valne for these gaps is shown in Figure ^as a function of 
(5/i with the NJL model resnlts from the previous section inclnded for com- 
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Figure 4: Comparison of the gaps as a function of djj, for the NJL model and 
single gluon exchange. 


parison. The results are also shown at 3 other hxed values of p in Figures 
I, I and 1^. These results show that both gaps consistently decrease with 
increasing over the whole range of momenta and the decrease is basically 
linear in 6p with a small quadratic component. 

If one takes the objects in square brackets from the equations for A^^{ud) 
and A^^{us) and in the NJL model, substitutes into the perturbative 
equation and examines the region away from the peak at / = 0 (/ >> Agg 
but not asymptotic values of 1), these terms are: 


A33M)(p) : 


Agg(Ms)(p) : 


5 

ICO 

ICO 

<] 

1 

18 

18 P 

5 

5 Agg J/i 

18 

36 P 


+ 

+ 0{5p^) 


(64) 

(65) 


Thus the effect of this region of the integration is to cause both gaps to 
decrease linearly in which is seen above. 

The other effect that one observes from Figure ^ is that the changes in the 
gaps in the NJL model are much more drastic than in perturbation theory. 
This is likely a result of the fact that, with the cut-off in momentum space 
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Figure 5: Gaps as a function of 6n in perturbation theory at p = 0. 

in the NJL model, the region of the integration near the Fermi momentum 
forms a much larger part of the complete integral, and so the relative effect is 
larger. In perturbation theory where one integrates over a larger momentum 
region within which different effects play a role, the relative change in the 
integral is much smaller. 

Summarizing the results, at /i = 500 MeV the dependence of A^^{ud) and 
on Sfi in first order perturbation theory is much weaker than in the 
NJL interaction model. More importantly, the dependence of A^slud) on Sfi 
actually changes sign. 

The implications of this change in sign could be signihcant. As Sfi in¬ 
creases, flavor is still locked to color except that s flavor and b color are 
distinguished from the other flavors and colors. Above Sfi = 2(A33)5^=o the 
two flavor color superconducting phase is favored and the gap in this case is 
larger than the ud gap in the CFL phase. If there is no intervening phase 
between the three flavor superconducting phase and the two flavor color su¬ 
perconducting phase, the fact that the A^^{ud) gap decreases suggests that 
the phase transition from the three to two flavor phases is first order. This 
analysis should be extended all the way to the three to two flavor phase 
transition point in order to verify this. As well, at this density, the inclu¬ 
sion of higher orders of perturbation theory must be examined to determine 
how they affect the results. Finally this reasoning ignores the possibility of 
intervening phases such as a three flavor equivalent of crystalline color super- 
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5[MeV] 

Figure 6 : Gaps as a function of (5/i in perturbation theory at p = 550. 


conductivity, which has been predicted in the two flavor case^], or CFL with 
meson condensation predicted in |Q in an effective Lagrangian model. It 
would be very interesting to study both of these phenomenon in the present 
framework. 

The results presented in this section were obtained using hrst order per¬ 
turbation theory at a moderate density where it should be considered a model 
and may not be accurate. One could speculate that at much higher values 
of the chemical potential where perturbation theory is dehnitely valid and 
the solution for A 33 is more strongly peaked at / = 0 [|T3, the differences 


between the perturbative results and the NJL model may not be as large as 
here. However, these results indicate that the dependence on 6fi should be 
weaker in perturbation theory than in the NJL model. 


5 Conclusion 


In this paper the general poles of the quasiparticle propagator in the case 
where the s quark has a different chemical potential than the other two 
quarks has been presented. These are a generalization of the poles given in 
0 . These poles are valid in a different model than |T^ and the relationship 


between the two models is not clear. 

These poles were used in numerical solutions of the gap equations in an 
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Figure 7: Gaps as a function of (5/i in perturbation theory at p = 1500. 


NJL model to demonstrate that they produce results for the gaps which are 
qualitatively consistent with BQ. 

The poles were also used in a numerical solution of the gap equation for 
A 33 (p) using perturbative single gluon exchange and to obtain fairly reliable 
estimates of all the gaps. The results obtained for = 0 are qualitatively 
and quantitatively similar to those in |]T^, which can be seen as support for 


both results. The results of are obtained in the two flavor case while the 
results presented here are obtained in the three flavor case. 

Finally, numerical results for the gaps were found as a function of 6fi out 
to 6fi = AbMeV ~ (A 33 ) 5 ^=o- The main conclusion is that perturbation the¬ 
ory at p = 500MeV predicts a decrease in A^^{ud) as a function of 5p where 
the NJL model predicts an increase in A^^{ud). This observation suggests 
that the phase transition from three to two flavor color superconductivity is 
hrst order. The second conclusion is that the dependence of both A^^{ud) 
and A^^{us) on Jp is much weaker than in the NJL model for p = 500MeV. 
Therefore, there are both qualitative and quantitative differences between the 
NJL model results and the leading perturbative results at moderate density. 
It is not clear which of these results should be closest to the real world. 

The weaker dependence of the gaps on Jp is likely to be true even at much 
higher densities, but how much different perturbative results might be from 
NJL results in this case remains to be seen. It is not clear if the increase of 
A^^{ud) with Jp will occur at higher p and this determination will be left for 
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a future paper. 

Analyzing the effect of shifting the chemical potential of the strange quark 
by 6fi is a hrst approximation to including the effect of a strange quark mass. 
This method has the beneht of separating out the difficulties introduced by 
fully including a strange quark mass in the gap equation. The knowledge 
gained in this research will be useful in an analysis where the s quark is 
given a non-zero mass at the level of the bare quark propagator which is the 
ultimate goal of this line of research]^. 

As well, the dependence of the gaps on Sfi presented in this paper may 
be interesting when compared with other papers that consider this effect. In 
it was shown in the two flavor case that up to 5/1 = Ao/a/2 , 5/r will have 
no affect on the gaps. In it is argued that in the three flavor case, 6fi will 


affect the values of the gaps and the number densities of the quarks but the 
densities of the different flavors of quarks will remain equal. This situation 
would mean that the superconducting ground state could remain electrically 
neutral without the addition of electrons which would drastically change the 
properties of this type of matter. This conclusion could presumably be tested 
using the results of this paper. 

This work can be extended in a number of ways. The analysis should be 
extended to = 2 (A 33 ) 5 ^=o corresponding to the phase transition to two 
flavors. The analysis should also be done at higher densities to determine 
how the results depend on p. Determine of higher order corrections to these 
results are necessary to determine the validity of the perturbative approach 
at this density. As well, the analysis of this paper could be extended is to 
non-zero temperature following [|^ . 
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